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Suppose several parties jointly possess a pure multipartite state, \ip). Using local operations 
on their respective systems and classical communication (i.e. LOCC) it may be possible for the 
parties to transform deterministically \ip) into another joint state \<fr). In the bipartite case, Nielsen 
majorization theorem gives the necessary and sufficient conditions for this process of entanglement 
transformation to be possible. In the multipartite case, such a deterministic local transformation is 
possible only if both the states in the same stochastic LOCC (SLOCC) class. Here we generalize 
Nielsen majorization theorem to the multipartite case, and find necessary and sufficient conditions 
for the existence of a local separable transformation between two multipartite states in the same 
SLOCC class. When such a deterministic conversion is not possible, we find an expression for the 
maximum probability to convert one state to another by local separable operations. In addition, we 
find necessary and sufficient conditions for the existence of a separable transformation that converts 
a multipartite pure state into one of a set of possible final states all in the same SLOCC class. 

Our results are expressed in terms of (1) the stabilizer group of the state representing the SLOCC 
orbit, and (2) the associate density matrices (ADMs) of the two multipartite states. The ADMs play 
a similar role to that of the reduced density matrices, when considering local transformations that 
involves pure bipartite states. We show in particular that the requirement that one ADM majorize 
another is a necessary condition but in general far from being also sufficient as it happens in the 
bipartite case. In most of the results the twirling operation with respect to the stabilizer group 
(of the representative state in the SLOCC orbit) plays an important role that provides a deep link 
between entanglement theory and the resource theory of reference frames. 

PACS numbers: 03.67.Mn, 03.67.Hk, 03.65.Ud 



I. INTRODUCTION 

Every restriction on quantum operations defines a re- 
source theory, determining how quantum states that can- 
not be prepared under the restriction may be manipu- 
lated and used to circumvent the restriction. Much of 
quantum information theory can be viewed as a the- 
ory of interconversions among such resources [T]. When 
several spatially separated parties sharing a composite 
quantum system, they are naturally restricted to act lo- 
cally on the system. Such a restriction to local opera- 
tions assisted by classical communication (LOCC) leads 
to the theory of entanglement, in which multipartite en- 
tangled states act as resources for overcoming the LOCC 
restriction. Hence, exotic multipartite entangled states 
play an important role in a variety of quantum informa- 
tion processing tasks. These include conventional [5] and 
measurement-based quantum computation [3], quantum 
error correction schemes [2], quantum secret sharing [4], 
quantum simulations [5], and in principle in the descrip- 
tion of every composite system consisting of more than 
one subsystem [6l [7] . 

The amount of information needed to describe n- 
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party quantum system grows exponentially with n, which 
makes it very difficult and almost impossible to clas- 
sify multipartite states by the degree of their resource- 
fulness (i.e. by their degree to circumvent the limita- 
tion of LOCC). The main reason is that the strength of 
the LOCC restriction increases with n, and, in general, 
the extent to which quantum resources (like multipar- 
tite entangled states) can be interconverted decreases as 
the strength of the restriction increases [8]. It can be 
shown HO] that for n > 3, the LOCC restriction leads 
to an uncountable number of incquivalent resources (i.e. 
multipartite entangled states). 

A conventional way to classify multipartite states is to 
divide the Hilbert space into sets of states that are all 
related to each other by invertible stochastic local oper- 
ations and classical communication (SLOCC) [9 1. More 
precisely, we say that tp and <fi belong to the same (in- 
vertible) SLOCC class if both transformations tp — > <f> 
and </>—>■ ip can be achieved by means of LOCC with 
non-zero probability. Hence, in general, if ip and <f> be- 
long to two distinct SLOCC classes, they are incompa- 
rable and essentially correspond to two deferent types 
of resources. While in three qubits there are only 6 in- 
equivalent SLOCC classes (see [5]), there are uncount- 
able number of inequivalent SLOCC classes in four (or 
more) qubits 0HO]. Roughly speaking, this fact implies 
that relative to the dimension of the Hilbert space, only 
a small number of transformations can be achieved by 
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local means (since the states involved must belong to the 
same SLOCC class). Nevertheless, the very few trans- 
formations that can be achieved locally are extremely 
useful for QIP tasks, as in the examples of transforma- 
tions involving the cluster states [3], graph states [IT], 
code states [2], and other novel multipartite entangled 
states [B]. 

A fundamental problem in quantum information is, 
therefore, to characterize all possible deterministic LOCC 
transformations among states in the same SLOCC class. 
That is, if ip and 4> belong to the same SLOCC class, 
what are the necessary and sufficient conditions that the 
transformation ip — »• <p can be achieved by LOCC with 
100% success (i.e with probability one)? Despite its im- 
portance, this problem was solved completely only for 
the case of pure bipartite states [6j[7]. In fact, even in 
the case of 3 qubits, very little is known about the char- 
acterization of local transformations among states in the 
same GHZ class [IS]. For higher number of qubits, or 
in higher dimensions, much less is known and there are 
only partial results for very specific SLOCC classes of 
states [TBIfTI]. 

One of the several difficulties encountered in the effort 
to solve this problem, is the non-elegant mathematical 
description of (the operationally motivated) LOCC pro- 
tocols in composite quantum systems consisting of more 
than two parties. This difficulty has been avoided in the 
literature by considering a larger set of transformations, 
such as the positive partial transpose preserving opera- 
tions |16j , or the non-entangling operations in |17) . How- 
ever, these sets of transformations are much larger than 
what can be achieved by LOCC, and therefore include 
many transformations that can not be achieved locally. 
A much smaller set, that on one hand is local and in- 
clude all LOCC transformations while on the other hand 
is mathematically simple to characterize, is the set of lo- 
cal separable transformations (SEP). 

SEP can be characterized as a set of trace preserving 
quantum channels such that for each channel in the set, 
there exists an operator sum representation (or Kraus 
representation) with the property that each Kraus ma- 
trix is a tensor product of n matrices acting on the n 
subsystems consisting the composite system (see the sec- 
tion [n] for the mathematical definition of SEP). In [TS] 
it was shown (somewhat counter intuitively) that SEP is 
strictly larger than LOCC, while in [TS] it was shown that 
for bipartite pure states we have effectively SEP=LOCC; 
that it, any transformation between two bipartite pure 
states that can be achieved by SEP can also be achieved 
by LOCC and vice versa. However, for mixed bipartite 
states this is no longer true |20) . 

In this paper we characterize the set of all transfor- 
mations taking one multipartite pure state to another by 
SEP. Clearly, unless the two states belong to the same 
SLOCC class, such a transformation is impossible. For 
each SLOCC class we show that there is a natural state 
(known also as a normal form [3_I]) that can be chosen to 
represent the class, and use the stabilizer group of this 



state to present our results. Unlike the stabilizer for- 
malism in quantum error correction schemes [], here the 
stabilizer group is not a subgroup of the Pauli group, but, 
in general, it is a subgroup of a much larger group, the 
group of all invertible matrices. 

When it is impossible to convert by SEP one state to 
another, we find an expression for the maximum prob- 
ability for the conversion. Our expression is in general 
hard to calculate and we therefore provide simple lower 
and upper bounds on the maximum probability of con- 
version. In addition, we generalize our results to include 
necessary and sufficient conditions for the existence of 
a separable transformation that converts a multipartite 
pure state into one of a set of possible final states all in 
the same SLOCC class. 

The paper is organized as follows. In section [TTJ we 
discuss some preliminary concepts and theorems about 
SLOCC classes, separable operations, and normal forms. 
In section |III| we define the stabilizer group and prove 
several of its properties. We also define the stabilizer 
twirling operation, and discuss the relation between the 
stabilizer group and SL-invariant polynomials. In sec- 
tion |IV| we define the associate density matrix and discuss 
its relation to the reduced density matrix. In Section [V] 
we present and prove our main results. Section VI is 



devoted to examples demonstrating how to apply our re- 
sults to specific cases. Finally, in section [VH| we end with 
conclusions. 



II. PRELIMINARIES 

The Hilbert space of n qudits is denoted by: 
H n = C dl <g> C d2 <g> ■ ■ • ® C d " . 
The SLOCC groups we will consider are the subgroups 



G = SL(di,C) 
G = GL(di,C) 



SL(ci 2 ,C)(§ 
)GL(d 2 ,C) 



8>SL(d„,C) 
®GL(d„,C) 



in GL(H n ). Here GL(d, C) is the group consisting of all 
dxd complex invertible matrices and SL(d, C) is a subset 
of these matrices with determinant 1. 

Viewing the SLOCC classes as orbits over the action 
of G, we define the orbit of a state ip £ H n to consists of 
unit vectors: 



Definition 1. The orbit of a state 
the action of G is defined by: 



£ T-L n under 



0,l 



urn 



g£G 



Note that the orbit above if unchanged if one replace 
G by G since we ignore global phase. Note also that 
any state in the orbit can be chosen to represent the 
orbit; that is, if ip' £ then 0^> — £)</,. Sometimes in 
literature the term "SLOCC class" is referred to a similar 
set as above, but slightly different, allowing the matrix g 
in the definition above not to be invertible. 
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A. Separable Matrices and Separable Operations 

The set of all states that can be prepared by LOCC 
is called the set of separable states. Here we define it 
without the normalization requirement of trace one. 

Definition 2. A matrix (or operator) p : H n —> W n is 
said to be a separable matrix or a separable operator if 
it can be written as: 

where all : C di C di are positive semi-definite ma- 
trices (operators). 

Later on, we will use the fact that any tensor product 
of matrices can be completed to the identity matrix by 
an addition of a separable matrix: 

Lemma 1. Let M — A\ % A-i ® • • • ® A n , where are 
dk x dk positive semi-definite matrices. Denote by A max 
the largest eigenvalue of M. Then, for p G M, the matrix 
I — pM is separable if and only if p < 1/A max . 

Proof. Clearly pM < I so I~pM is positive semi-definite. 
To see that it is also separable denote Xj to be the max- 
imum eigenvalue of Aj, denote Ij to be the identity 
matrix on subsytem j, and w.l.o.g take p = 1/A max = 
(A1A2 ■ • • A„) . Then, the matrix Ij — j-Aj is positive 
semi definite, and therefore Ij is a sum of two positive 
semi definite matrices Ij = (Ij — j-Aj) + j-Aj. Substitut- 
ing this into I = I\ ® I2 <8> • • • I n completes the proof. □ 

We now give the precise definition of the set of opera- 
tions that will be used in this paper. 

Definition 3. Let B(H n ) be the set of density matrices 
acting on H n , and let A : B(H n ) — > B(H n ) be a com- 
pletely positive trace-preserving map. Then, we say that 
A is a separable super-operator if there exists an operator 
sum representation of A such that for any density matrix 
peB(n n ), 

A(p) = Y,M kP Mt, (1) 

and 

M k = 4 fe) ® ® • • • ® ^l fe) » 

where are di x d{ complex matrices. 

We denote the set of all separable superoperators by 
SEP. In the particular case, when considering only pure 
bipartite states, LOCC is effectively the same as SEP, 
but in general, LOCC is a strict subset of SEP. 



B. The Critical Set 

Let g be the Lie algebra of G contained in End(H n ). 
We define the set of critical states as 

Crit(Hn) = {4> G H n \ (0|A|0) =0,l£ 0}. 

The GHZ states [9], cluster states [3], graph states [IT], 
and code states [2] (such as the 5-qubit codes state, the 
Steane 7-qubit code state, and the Shore 9-qubit code 
states) are all critical. In the bipartite case, with H2 = 
C d ®C d , the maximally entangled state is the only critical 
state up to a local unitary matrix. As we will see shortly, 
almost all states belong to an orbit of some critical state. 

The Kempf-Ness theorem 22 in the context of this 
paper says (the only hard part of the theorem is the "if 
part of 3): 

Theorem 2. Let 4> G U n and let K = SU{d 1 )®SU{d 2 )® 
■■■(g) SU(d n ), then 

1. 4> G Crit(H n ) if and only if \\g4>\\ > \\4>\\ f or a ^ 

9 e G. 

2. If $ e Crit(Hn) and j € G then \\g<j>\\ > \\(/)\\ 
with equality if and only if g<p G K<p. Moreover, if g is 
positive definite then the equality condition holds if and 
only if g<p = <f>. 

3. If <f> G H n then Gift is closed in H n if and only if 

Gcj)r\Crit(n n ) f 0. 

The dimension of the set Crit(H n ) in the case of n 
qubits (i.e. d\ = di = ... = d n = 2) is 1 for n = 3 (up 
to a local unitary matrix, only the GHZ state belongs to 
Crit{U n )) and 2" - 3n for n > 3. 

The following theorem says that critical states are 
maximally entangled. 

Theorem 3. US WM Let i> e Hn- Then , i> & Crit(U n ) 
if and only if the local density matrices are all propor- 
tional to the identity (i.e. each qudit is maximally entan- 
gled with the rest of the system). 

If an SLOCC orbit contains a critical state i/> c , then 
ij} c can be regarded as the maximally entangled state of 
the orbit O^. In such cases, we will take ip c to represent 
the orbit and denote the orbit by 0^ c . The corollary 
below implies that almost all SLOCC classes can be rep- 
resented by a critical state. 

Corollary 4. \2S^ The set of all states of the form 
gip/\\gip\\, where g £ G and ip G Crit(J-L n ), is dense in 
H n . 

Moreover, part 2 of the Kempf-Ness theorem above 
implies the uniqueness of critical states; that is, up to 
local unitary matrix there can be at most one critical 
state in any SLOCC class O^. That is, critical states 
are the natural representatives of SLOCC orbits, and are 
the unique maximally entangled states in their SLOCC 
orbits. 
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The main results of this paper are expressed in terms of 
the stabilizer group of the state representing the SLOCC 
orbit, and the multipartite generalization of the reduced 
density matrices, which we call the associate density ma- 
trices (ADMs). In the next two sections we define and 
prove useful theorems about stabilizer groups and ADMs. 
These will be fruitful when we present our main results. 



III. THE STABILIZER GROUP 

In section [V] we will see that the stabilizer group natu- 
rally arrises in the analysis of separable operations. How- 
ever, unlike the stabilizer formalism of quantum error 
correcting codes [2], in general, the stabilizer group con- 
sidered in this paper is not a subgroup of the Pauli group. 
Instead, it is a subgroup of G, which is a much bigger 
group than the Pauli group. 

Definition 4. The stabilizer group of ip £ % n is a sub- 
group of G defined by 

Stab(^) = {g £ G\ gi, = V^} . 

Note that for g £ G we have Stab(g?/;) = gStab^)^ -1 . 
That is, the stabilizer groups of two states in the same 
SLOCC class are related to each other by a conjugation. 

Proposition 5. // the stabilizer group Stab(ip), for ip £ 
T-L n , is finite, then there exists a state </> £ O^p such that 
Stab(4>) C U{di) <S> U(d 2 ) ® • • ■ (8> U(d n ). 

Proof. The proof is by construction. Denote by 
Stab(ijj) = {Sk}k=i,2,...,m, where m is the number of el- 
ements in Stab(%jj). Since Stab(^p) is a subgroup of G, it 
follows that the elements S k can be written as: 



S k = A 



(i) 



.4 



(2) 



.(«) 



where £ GL(g?;,C) and I = 1,2, ...,n. Clearly, for 

a fixed I, the set {^4^}/=i,2,...,n form a subgroup of 
GL(d;,C). This subgroup is equivalent (up to a con- 
jugation) to a unitary group. To see that, denote, 



1/2 



\k=l j 



(I) - 



Clearly, £ GL(d ; ,C) and the elements U, 

A^A^A^) -1 form a unitary representation of the sub- 
group {A^}. Now, denote by 



A (1) <g> A (2) ® • • • <g> A (n) 



<8> U, 



(n) 



(2) 



u k = ujp ® ui 2) c 

By construction we have 

U k A\i/j) = A|V> ■ 

Hence, the stabilizer group of (j> = A\tji) /\\A\tjj)\\ consists 
of unitaries. □ 



Note that for generic states of 4 or more subsys- 
tems [23] the stabilizer group is finite and therefore the 
proposition above can be applied to most (i.e. a dense 
subset) of the states in W n . We now show that if the or- 
bit of a state with finite stabilizer contains a critical state 
then the stabilizer of the critical state is unitary. Since 
most orbits contains a critical state, it is most natural to 
choose the critical states to represent the SLOCC orbit 
(see Fig.l). 

Proposition 6. If the stabilizer group Stab(ip) for ip £ 
Jin is finite, and there exists a critical state ?p c in the 
orbit O^, then Stab(tp c ) C U(di) ® U{d 2 ) ® • • • ® U(d n ). 

Proof. Let g £ Stab(tp c ) C G; i.e. g\ip c ) — IV'c) • In its 
polar decomposition g = up, where u £ U(d\) ® U{d2) <g> 
• • • ® U (d n ) and p £ G is a positive definite matrix. Note 
that ||V'c|| — llsV'cll = IbV'cll- From theorem [2] (Kempf- 
Ness), \if> c ) being critical implies that UpV'cII ^TlV'cll with 
equality if and only if p\ip c ) = \ip c ) ■ This implies that p £ 
Stab(ip c ). Moreover, since tp c £ there exists h £ G 
such that Stab(ip c ) — hStab(i/))h~ 1 . From proposition [5] 
Stab(ip) is itself conjugate to a unitary subgroup. Thus, 
p is similar to a unitary matrix. However, p is positive 
definite matrix and since it is similar to a unitary matrix 
it must be an identity. □ 

Even though for most states in "H„ the stabilizer is 
a finite group, for many important states such as GHZ 
states, and other graph states, the stabilizer group is not 
finite. For such non-generic states we show that the sta- 
bilizer can not even be compact. 

Proposition 7. Let <fi £ H n be such that Stab(<fi) is 
compact then Stab(<f>) is finite. 

Proof. The stabilizer of an element is a linear algebraic 
group over C. The identity component of this group is 
therefore also a linear algebraic group over C. The only 
connected, linear algebraic group that is compact is the 
group with one element. □ 

When the stabilizer group of a state ip £ H n is not 
finite and therefore non-compact, it would be very useful 
to work with a compact subgroup T C Stab(xp), defined 
as the intersection of Stab(ip) with U(di) <E> U (d 2 ) ® • • • ® 
U(d n ). We will discuss it in more details in Sec VD and 
Sec ED 



A. The Stabilizer Twirling Operation 

In this subsection we define the stabilizer twirling op- 
eration that will be used quite often when we discuss our 
main results in the next section. 

Definition 5. Let ip £ W n be a state with a finite uni- 
tary stabilizer. The ^ta^i/^-twirling operation, G(-), is 
defined by 



(3) 



k=l 
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for any positive semi-definite operator a acting on % n . 

In Ref. [23] it was shown that the twirling operation 
Q can be factorized in terms of the irreps of the group 
(in our case Stab(ip)). In the following we make use of 
techniques introduced in Ref. |24j , to which the reader is 
referred for more details. Under the action of the uni- 
tary representation of a compact group Stab(ip), a finite 
dimensional Hilbcrt space factorizes as follows 



where q labels the irreps of Stab(ip), M. q is the qth rep- 
resentations space, and N q is the gth muliplicity space. 
The Stab(ip)-twix\ix\% operation has the form 

g (a) = V v Mq ® id^ 9 (n w ,rfi Wt ) , (4) 



where Tln q is the projector onto T-L q , idjv is the identity 
map on J\f q , and T>M q is the completely decohering map 
on M. q ; i.e. it denotes the trace-preserving operation that 
takes every operator on the Hilbert space Ai q to a con- 



stant times the identity. In section VA we will see that 



this factorization enables the elegant classification of all 
states <f> to which ip can be transformed dctcrministically 
by separable operations. 

As we pointed out above, if the stabilizer group is not 
finite, then it is also non-compact and therefore there is 
no Stab(ip)-twirYnig. However, for the compact subgroup, 
T C Stab(ip) as defined above, there is an analog for 
the twirling operation. We extend the definition of the 
stabilizer twirling operation (see Eq.([3])) to the compact 
group T. 

Definition 6. Let ip £ T~L and assume Stab(ip) is non- 
compact. Let T be the intersection of Stab(ip) with 
U(d\) ® U(d,2) <E> • • • <E> U(d n ). Then, for any positive 
semi-definite operator a acting on T-L n , the T- twirling op- 
eration, T(c), is defined by 



T(ct) = J dt $ot 



(5) 



where t £ T and dt is the Haar measure over T. 



The factorization form in Eq.Q holds for the T- 
twirling as well. This T- twirling operation will become 
useful whenever the stabilizer group is non-compact. The 
simplest example of that is the 3 qubits GHZ state (see 
secfVTl. 



B. SL- Invariant Polynomials 

In this subsection we discuss the quantification of en- 
tanglement in terms of SL-invariant polynomials and dis- 
cuss their relationship to the stabilizer group. 

An SL-invariant polynomial, f(tp), is a polynomial in 
the components of the vector ip £ H n , which is invariant 



under the action of the group G. That is, f(gip) = f(ip) 
for all g £ G. In the case of two qubits there exists only 
one unique S'L-invariant polynomial. It is homogeneous 
of degree 2 and is given by the bilinear form (ip, ip): 



f 2 (ip) = (ip,ip) = (r 



ip eC 2 ®C 2 



Its absolute value is the celebrated concurrence [23] , 

Also in three qubits there exists a unique S'L-invariant 
polynomial. It is homogeneous of degree 4 and is given 
by 



U(iP) = dot 



OPQj^o) {<P0,<Pl) 

(PuVo) (<Pi,<Pi) 



where the two qubits states ifi for i = 0, 1 are defined 
by the decomposition \ip) = |0)|</?o) + an d the 

bilinear form (tpi,ifj) is defined above for two qubits. 
The absolute value of f± is the celebrated 3-tangle . 

In four qubits, however, there are many SL-invariant 
polynomials and it is possible to show that they are gen- 
erated by 4 SL-invariant polynomials (see e.g. [TU] for 
more details and refrences). 

Any SL-invariant polynomial is generated by a fi- 
nite number of homogeneous SL-invariant polynomials. 
Moreover, for some states, all SL-invariant polynomials 
vanish. For example, the W state 

| WO = (|100) + |010) + |001))/V3 

and (its generalization to n-qubits) has the property that 
g t \W) = t\W), where g t = diag{t, t' 1 }^ 3 and ( £ I is 
non-zero. Note that g t G SL(2,C)® 3 . Thus, the value of 
any homogeneous SL-invariant polynomial, /, is zero on 
\W) (and its generalization to n qubits) since 

f{\W)) = f{g t \W)) = f(t\W))=t k f(\W)), 

where k £ N is the degree of /. Nevertheless, the set of 
states for which all SL-invariant polynomials vanish is of 
measure zero and is called here the null cone. Since the 
null cone is of measure zero, for most states there exists 
at least one homogeneous SL-invariant polynomial that 
does not vanish on them. The following proposition will 
be useful for such states. 

Proposition 8. Let ip £ Ti. n and assume ip is not in the 
null cone. Let k £ N be the degree of some homogeneous 
SL-invariant polynomial that is not zero on ip. Then, 



Stab(iP) C G k 



where 



Gk 



[g' £ G\g' = e^g ; 1 < m < k ; g £ 67 } 



is a subgroup of G. 

Proof. Let / be an homogeneous SL-invariant polynomial 
of degree k such that f(ip) ^ 0. Let g' be an element in 
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Stab(ip). Since Stab(ip) C G, there exists a s C such 
that g' = ag, where g G G. Thus, 

fW = /(ffV) = / (aflVO = «"/ (#) = a"/ WO , 

where the first equality follows from the fact that g' € 
Stab(ip), whereas the last equality follows from the fact 
that g € G and / is SX-invariant polynomial. Since 
f(i/j) 7^ it follows that a fe = 1. This completes the 
proof. □ 

Remark. If k, in the proposition above, is equal to the 
dimension, rfj, of the subsystem I (here Z = 1, 2, n), 
then G k = G. This is because e t27rm / d, g € G for all 
g E G, I = 1,2, ...,n and m = 1,2, In particular, 

for systems of n qubits G 2 = G. 

In the following corollary, we show that for most states 
in H n (i.e. a dense subset of the stabilizer group is a 
subgroup of G (which is a smaller and somewhat simpler 
group to work with than G). 

Corollary 9. Let fk and fk> be two homogenous SL- 
invariant polynomials of degrees k and k 1 , respectively. 
Let ip € H n , and assume that fk{ip) ^ and fk>(ip) ^ 0. 
Denote by I = (k,k') the greatest common divisor (gcd) 
of k and k' . Then 

Stab(ip) C Gi . 

This corollary follows directly from the proposition 
above. Note also that since most states in T-L n has many 
non-vanishing SL-invariant homogenous polynomials, it 
follows from the corollary above that for these states the 
stabilizer group is a subgroup of G. For example, consider 
the case of n-qubits with even n. For such systems their 
exists a homogeneous SX-invariant polynomial of degree 
2. It is possible to show that for a dense set of states in 
T-L n this polynomial has non-vanishing value. Thus, for 
this dense set of states the stabilizer is a subgroup of G2 
which is equal to G for qubits. 

Despite the above corollary, in many important cases 
(such as the 3 qubits GHZ state) the stabilizer is not a 
subgroup of G. For these cases the following corollary 
will be useful. 

Corollary 10. Let ip g H n and define St(ip) to be the 
intersection of Stab(ip) with G; that is, 

Sty) = {g e G\ gi, = ^} . 

If there exists an homogeneous SL-invariant polynomial 
of degree k that is not vanishing on ip, then the quotient 
group Stab(ip) / St(ip) is a cyclic group of order at most 
k. 




SLOCC CLASSES 



FIG. 1: The state space H n can be described heuristically 
as a horizontal cylinder. The horizontal dimension represents 
the uncountable number of inequivalent SLOCC classes. Each 
SLOCC orbit of states is described by a circle perpendicular 
to the horizontal line. The states (pi and <fi2 are the repre- 
sentatives of two inequivalent SLOCC orbits. The action of 
g £ G on ipi and ip 2 generates <p'i and tp 2 , respectively. There- 
fore, the ADMs of both ip\ and (p' 2 are proportional to g^g. 
From the results in section lV Al it follows that if the stabilizers 
of cpi and ip 2 are equal, then the transformation ipi — > Lp[ is 
possible by SEP iff ip 2 — » V?2 is possible by SEP. 



IV. THE ASSOCIATE DENSITY MATRIX 

In this section we associate density matrices (not nec- 
essarily normalized) with states in an orbit O^. We 
call such density matrices associate density matrices 
(ADMs). We will see that the ADM plays a similar role 
to that of the reduced density matrices, when considering 
local transformations that involves pure bipartite states. 

Definition 7. Let ip € H n and let cf> € be another 
state in the SLOCC orbit of ip. Thus, there exists g 6 G 
such that 4> = gij)/\\gip\\. An associate density matrix, 
p^{4>), of 4> with respect to ip, is defined by 

Note that the associate density matrix is positive def- 
inite and that for s € Stab(V'), both p^(<f>) and s' p^,{(p)s 
are associate density matrices of (p with respect to ip. 
Hence, p^((p) is defined up to a conjugation by the stabi- 
lizer group of ip. Furthermore, note that p^,(</>) = / if and 
only if ip is related to (p by a local unitary transformation. 

Proposition 11. Let ip,<p E H n be two states in the 
same orbit (i.e. <p = gip/\\gip\\, for some g G G), 
where Stab(ip) C U(di) <g> U{d 2 ) ® ■■■ ® U(d n ). Then, 
Stab((p) C U(di) U(d 2 ) (8 • • • U{d n ) if and only if 

[p-)p(4>), u] = for all u £ Stab(ip) . 

Proof. Since <p — 9^1 1 1 9i> 1 1 we have Stab(<p) = 
gStab{ip)g~ 1 . That is, if s e Stab{<p) then there exists 
u e Stab(ip) such that s = gug~ x . Thus, s is unitary if 
and only if 

s f s = (s 1 )" 1 " 1 S 1 9 U 9^ = 1 u^g^gu = g^g . 



This corollary is a direct consequence of proposition [8] 



This completes the proof. 



□ 
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The relation with the reduced density matrix 

Consider the Hilbert space C d ®C d , and the maximally 
entangled state (it is unnormalized for simplicity) 

d-l 



where {\i)} is an orthonormal basis in C d . Clearly, in 
this case the orbit consists of all states in C d <X> C d 
with Schmidt number d. 



The stabilizer group of ip is given by (see Sec. VIA for 
the proof): 



abty) = {s- 1 ®S T S e GL(C d )} , 



St 



where S T is the transpose matrix of S. 
Any state <j> € can be written as 



(6) 



A 2 \ip}\ 



A 1 ®A 2 



where A\®A 2 £ G and in this subsection G — SL(d, C)(g> 
SL(d,C). However, due to the relatively big stabilizer 
group, we can always assume without loss of generality 
that A 2 = I. Thus, the associate density matrix (which 
is defined up to conjugation by an element from the sta- 
bilizer group) can be written as: 



Am 



\Ai® I\tp)\\ 



where p r = Tr b|<^) '{4>\ IS the reduced density matrix. 
Thus, in the bipartite case, the associate density matrix 
is essentially equivalent to the reduced density matrix. 
As we will see in the next section, in the multipartite 
case, the associate density matrix plays the same role 
as the reduced density matrix plays in bipartite LOCC 
transformations. 



V. MAIN RESULTS 

In this section we present our main results. As we have 
seen in the previous sections, for a dense set of states in 
Wn the stabilizer group is finite. Hence, in the first three 
subsections we assume explicitly that the stabilizer is fi- 
nite, while our results for infinite stabilizer groups are 
presented in the last subsection. More specifically, in 
subsection |V A| we generalize Nielsen majorization thc- 

we 



orem to multi-partite states, and in subsection VB 



introduce necessary and sufficient conditions for non- 
deterministic SEP transformations. Subsection IV CI is 
devoted to the calculation of the maximum probability 
to convert one state to another in the case where deter- 
ministic SEP is not possible. Finally, in section [VD| we 
generalize our theorems for the case in which the stabi- 
lizer group is non-compact. 



A. Generalization of Nielsen Majorization 
Theorem to Multipartite states 

In [27] it was shown that a transformation | -01 } AB — » 
\ip 2 ) AB , between two bipartite states in C™ <g)C m , is pos- 
sible by LOCC if and only if there exists a set of prob- 
abilities {pk} and a set of unitary matrices {Uk} such 
that 



ft 



fc 



PkUipfUk , 



where pf = TrslV'i) (V^l — 1) 2) are the reduced den- 
sity matrices. Moreover, it was argued that the relation 
in the equation above hold true if and only if the eigen- 
values of pi are majorized by the eigenvalues of p 2 ; i.e. 
Pi Pi - We now generalize this theorem to multipartite 
states. 

Theorem 12. Let |-0i> = S , i|^)/||s , iV>ll and \ip 2 ) = 
<72|V')/II<?2' ! / ; II be two entangled states in the same orbit 
O^p, where the orbit representative \ip) G % n is chosen 
such that Stab(\ip)) is a finite unitary group (see propo- 
sition^. Let m be the cardinality of Stab(ip). Then, 
\tpi) — > 1^2} by SEP if and only if there exists a set of 
probabilities {pk} such that 



^2pkUlp 2 U k = pi , 



(7) 



fc=i 



where p t = p^{ipi) for i = 1,2 and Uk € Stab(ip). 



Remark. The condition on p\ and p 2 in the Eq.([7| is 
the multipartite generalization of the condition obtained 
in |27j for the bipartite case. In particular, Eq.([7]) implies 
that pi -< p 2 , i.e. the eigenvalues of p\ are majorized 
by the eigenvalues of p 2 . In the bipartite case, this ma- 
jorization condition is both necessary and sufficient. The 
reason for this is that in the bipartite case the stabilizer 
group is big and contains all unitary matrices, whereas 
in the multipartite state the stabilizer usually consists of 
only a finite number of unitary matrices Uk- Hence, in 
the multipartite state the majorization condition is far 
from being sufficient. 

Proof. Consider the separable operation given in Eq.Q. 
If (ipi\) = \ip%) (ip2 1 then there must exist complex 

coefficients such that 

M k \ipi) = a k \^ 2 ) 

for all k. Since the separable operation A is invariant un- 
der the multiplication of the operators by phases, we 
can assume without loss of generality that are real and 
non-negative. We therefore denote au = \Jp~k~, where pk 
are non-negative real numbers. From the above equation 
and the completeness relation J^k -^fc-^fc = ^ it follows 
that YskPk = 1- 



Now, since \ip\) 
92\i>)/\\g2ip\\ we get 

IIMIL 
HfliV'llv^fe 



SilV'Vllfl'iV'll and \tp 2 



9z l M k gM) = m 



Hence, 



llfla^ll 
\gii>\\y/Pk 



g 2 1 M k g 1 = U k 



where U k € Stab(tp). In the above equation we used 
the fact that tp is n-way entangled (in the sense that the 
determinant of each of the n reduced density matrices ob- 
tained by tracing out n—1 qudits, is non zero); otherwise, 
M k may have zero determinant. The above equation may 
be rewritten as: 



\92ip\\ 



/pk92U k g 1 



(8) 



The completeness relation ^2 k M k Mk — I gives: 

m 

^2,PkUlp 2 U k = px , 

k=l 

where we substitute pi = glgi/\\giip\\ 2 for i = 1, 2. This 
provide the first direction of the proof. The other direc- 
tion follows the exact same lines backwards. □ 

Taking on both sides of Eq.([7]) the S , ta6('i/ , )-twirling 
operation gives the following simple necessary condition 
which we summarize in the following corollary: 

Corollary 13. Notations as in theorem\T^ If the trans- 
formation \ipi) — > \ip2) is achievable by SEP then 



G(p 2 ) = G(pi) 



(9) 



Moreover, if p\ is invariant under the action of the 
Stab(ip) -twirling operation, that is, if 



GiPi 



Pi 



then the condition in Eq. |$J) is both necessary and suffi- 
cient. 

The corollary above implies in particular that the 
transformation ip — > -02 is achievable by SEP if and only 
if 



G(p 2 



(10) 



We have seen in the previous sections that if i\) G H n is 
both critical and generic, then Stab(ip) — {U k }k=i 2,...,m 
is a finite unitary group. Hence, from Eq. ( 10 ) and 
Eq. @ we can characterize all the states <f> that can be 
obtained from a critical generic state ip. We summarize 
this characterization in the following corollary: 

Corollary 14. Let ip € T-L n be a state with a unitary 
stabilizer group, and let <p be a state in the orbit O^p. 
Then, the transformation ip — > <p can be achieved by SEP 
if and only if 

Tr Mq [n« 9 p^( ( / , ) n «,] = id A/-, V q . 



B. Non Deterministic Transformations 

Transformations involving quantum measurements are 
usually not deterministic. In this subsection we provide 
necessary and sufficient conditions for the existence of a 
local procedure that converts a finite multipartite pure 
state into one of a set of possible final states. For the 
bipartite case, this problem was solved in [25- In partic- 
ular, a minimal set of necessary and sufficient conditions 
where found for the local conversion of a bipartite state 
■ 1B to the ensemble of bipartite states {pj, \4>j) AB }', 



i.e., \ip) AB is converted to 
The conditions are 



\AB 



with probability pj. 



AB\ 



where {E k } are the Vidal's entanglement monotones |29j : 



AB\ 



for k = 1, 2, n 



(11) 



where Ai > A 2 > • • • > A„ are the ordered eigenvalues of 
the reduced density matrix p A = Tr b\iP) ab (ip\- The fol- 
lowing theorem generalizes these conditions to the mul- 
tipartite case. 

Theorem 15. Let <p G Ti. n be a normalized state in the 
orbit 0^, where the orbit representative \ip) is chosen 
such that Stab(ip) consists of m unitary matrices (see 
proposition^^. Let {pj,<j>j} be an ensemble of states also 
in the orbitU^. Then, the transformation cf> — > {pj,<fij} 
is achievable by SEP (i.e. there exists a separable oper- 
ation taking ip to (pj with probability pj) if and only if 
there exist probabilities {pj k } such that: 



Eft* 
fe=i 



Pj 



3 k=l 

where p = p^{4>j), Pj = p^(<j)j), and U k € Stab(ip). 



(12) 



Remark. The relation between p and pj in Eq.(12| is 
the multipartite generalization of the condition obtained 
in [55] for the bipartite case. In particular, Eq.(12| im- 
plies that 



£ fe (|0»>]>>^(l^» > 



(13) 



for all k = 1, 2, dim"H n . On multipartite states the 
functions {E k } are defined by 



dim«„ 



E k (|0)) = A ™ ' for k = 2 > ->dimH„ (14) 



m—k 



where Ai > A 2 > ■ ■ • > XdixaU n are the ordered eigen- 
values of the ADM of \<fi). In the bipartite case, these 
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conditions are both necessary and sufficient since the sta- 
bilizer group contains all unitary matrices, whereas in the 
multipartite case the stabilizer usually consists of only a 
finite number of unitary matrices Uk- Hence, in the mul- 
tipartite case the conditions in Eq.( 13 1 are far from being 
sufficient. 

Proof. The proof follows similar lines as the proof of thc- 
T2l Consider the SEP A = £\ Aj , where 



K 3 { P ) = Y. M ^P M }k 
fe=i 

for any density matrix p acting on T-L n . If the transfor- 
mation <p — > {pj,4>j} is achievable by such a SEP, then 

Aj(MM)=Pil&)<&l ■ 

Thus, there exists complex coefficients cijk such that 

M jk \<t>) = ajklfa) 

for all j and k. Since the separable operation A is in- 
variant under the multiplication of the operators Mjk 
by phases, we can assume without loss of generality 
that cijk are real and non-negative. We therefore denote 
a.jk = y/Pjk, where pjk are non-negative real numbers. 
From the above equation it follows that J2kPjk — Pj- 

Now, since \<j>) = g\^) /\\gtp\\ and \(f>j) = 9j\ip) /\\gjip\\ 
for some g, gj £ G, we get 



\m\\y/Pjk 3 



Hence, 



/Pjk9i 1 M jk g = S k 



where Sk € Stab(ip). In the above equation we used 
the fact that ip is n-way entangled (in the sense that 
the determinant of each of the n reduced density matri- 
ces, obtained by tracing out n — 1 qudits, is non zero); 
otherwise Mjf. may have zero determinant. The above 
equation may be rewritten as: 

\\gij)\\ — _j 
Thus, the completeness relation V fe M^ k Mjk = I gives: 



J2J2 P 3 kU kP3 U h = P > 



k=l 



where we substitute p = Pip{4>j) and pj = g]gj /\\gj4>\\ 2 ■ 
This provide the first direction of the proof. The other 
direction follows the exact same lines backwards. □ 



Taking on both sides of Eq.(12) the 5ia6(^)-twirling 
operation gives the following simple necessary condition 
which we summarize in the following corollary: 



Corollary 16. Notations as in theorem\T5[ If the trans- 
formation \<f>) — > {pj, \ 4>j}} is achievable by SEP then 



^2pjG(Pj) = Q{p) ■ 



(15) 



Moreover, if p is invariant under the action of the 
Stab{%p) -twirling operation, that is, if 

G(p) = P, 

then the condition in Eq. ft is both necessary and suf- 
ficient. 

The corollary above implies in particular that the 
transformation ip — > {pj , <pj } is achievable by SEP if and 
only if 

^2pjG(pj) =1 ■ 



C. Probabilities 

From the theorems above we see that when the conver- 
sion of one state to another is possible with some prob- 
ability, i.e. the states in the same SLOCC class, this 
probability is usually less than one. A natural question 
to ask is what is the maximum possible probability, P max , 
to convert one state to another if we know that the two 



states are in the same SLOCC class. Eq. ( 13 1 provides an 
immediate upper bound on P max , which we summarize 
in the following corollary. 



Corollary 17. Lettpi and ^2 be as in theorem 12 Then, 



-Pmax W'i W < min <^ 



(16) 



where k = 1, 2, dim'Hn and E k are the functions de- 



fined in Eq.(H) 



Proof. The procedure converting \ipx) to \1p2j with prob- 
ability P max also converts IV'i) to some other states with 
the remaining probability 1 — P ma x- The value of E k 



for these states is non-zero and therefore from Eq.(13) it 
follows that P max < Ekiip^/Ek^)- 



□ 



In [35] it was shown that for the bipartite case the 
upper bound in Eq.( 16 ) is always achievable. In the mul- 



tipartite case, unless ip2 is critical, this is usually not the 
case and as we will see below in some cases it is even 
possible to find better upper bounds. In particular, it 
is possible that for all k, Ek(ipi) > -Eft (^2) and yet the 
probability P max is strictly less than 1 (or even arbitrarily 
close to zero). 

In the following we first derive a mathematical expres- 
sion for the maximum probability. Since in general the 
expression can be difficult to calculate (at least as hard 
as determining if a multipartite density matrix is sep- 
arable or not), we simplify it dramatically for the case 
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in which the states ipi or ip2 are critical. Then, for the 
more general case, we find lower and upper bounds for 
-Pmax- These upper and lower bounds remain the same if 
we replace SEP by LOCC. 



1. General formula for Pmax 

Theorem 18. Let ip, ipi, ip 2 , Pi 1P2 and m be as in 

For any set of probabilities {Pk}™=i with 
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theorem 

Y^k=iPk < 1; denote 



A Vl-^2 ({Pk}) = Pi 



k=l 



where Uk £ Stab(ip). Then, the maximum probability to 
convert tpi to tp 2 is given by the following expression: 

= max|y^p fc | A^^ 2 ({pk}) is separable^ . (17) 

That is, the maximum is taken over all sets of probabili- 
ties {pk} such that J2T=iPk — 1 an d A^ )1 ^ 2 ({pk}) is a 
positive semi-definite separable matrix. 

Proof. Any separable transformation that takes tpi to tp 2 
with maximum probability P ma x, takes tpi to some other 
states that are not tp 2 with the remaining probability 1 — 
Pmax- However, since any state can be transformed to the 
product state |00...0) with 100% success , we can assume, 
without loss of generality, that the separable map that 
takes tpi to tp 2 with maximum probability P ma x, takes tpi 
to the product state |00...0) with probability 1 — P max - 
A separable transformation with such property, can be 
written as A = Ai + A2, where 



fc=i 



is a separable map that takes -01 to tp 2 with probability 
Tr[Ai(Vi)] and 



Aa(-) = J2 M k(-)Ml 



is a separable map that takes tpi to 1 00. . .0) with proba- 
bility Tr [A 2 Oi)] = 1 - Tr [A^i)]. The form of M k for 
k < m is given in Eq. ^ . Note that Eq. ([8| also implies 
that there are at most m (i.e. the cardinality of stabi- 
lizer group) distinct Kraus operators in the (separable) 
operator sum representation of Ai. If there are strictly 
less than m Kraus operators in the representation of Ai , 
then we can complete it to m Kraus operators by adding 
zeros (i.e. taking some of the p k s in Eq.(|8]) to be zero). 
With these notations, the map Ai converts tpi to ip 2 with 
probability 



Tr[A 1 (^ 1 )]=£ 



Pk 



Now, the completeness relation for the tace-preserving 
CP map A, can be written as, 

m 

1 = E M t M * = E M l M * + E M i M * 



k=l 



fe>m 



m _ 



fc=i 



where we have used Eq. |8]). Hence, the probabilities 
{pk} for Ai are chosen such that 



({Pk}) 



1 



11511 



E M i M * 



91 



\k>m 



That is, ({pk}) is a positive semi-definite sepa- 

rable matrix. This completes the proof. □ 



The requirement in theorem 18 that A^-^j ({Pk}) is 
a separable matrix implies that Q(A^ )1 ^.^ 2 ({pk})) is also 
a separable matrix. This leads to the following necessary 
condition: 

Corollary 19. Notations as in theorem [7&[ If P m = 
Pmaxiip 4>) then the matrix 

Q{px) - P m Q (p 2 ) (18) 

is separable. 

As we will see shortly, this leeds to an upper bound on 

p 

2. A simple formula for states with unitary stabilizer 

If the state tp 2 in theorem 18 is the critical state tp, then 
p 2 = I and the expression A^ 1 ^ f ^ 2 ({pk}) in theorem 
is given simply by 



18 



A V>i^2 ({Pk}) = Pi 



k=l 



Pk I 



fc=l 



Now, from the similar arguments as in lemma[T]it follows 
that A i p 1 ^ 2 {{p k }) is separable if and only if X^fcLi-Pfc 
is not greater than the smallest eigenvalue of p\. We 
therefore have the following corollary: 

Corollary 20. Let tp € rl n be a critical state and let 
tpi G Oip. Then, the maximum possible probability to 
convert tpi to tp is given by 

-Pnax^l -> VO = A m in[Pl] , 

where A m j n [pi] is the minimum eigenvalue of p\. 

Note that in this case the upper bound given in corol- 
lary[l6]is saturated by this value. That is, the probability 
to convert a state to the maximally entangled state of the 
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SLOCC orbit (i.e. the critical state) is given by the same 
formula that it is given in the bipartite case. 

The expression in theorem [18] can also be simplified if 
0i is critical or if the stabilizer group of ip% is unitary; 
i.e. if ipi is also a natural representative of the orbit O^p. 
Since most states are generic, their stabilizer group is 
finite, and if they are also critical then from proposition]!}] 
it follows that their stabilizer group is unitary. For such 
states the following corollary applies. 



Corollary 21. Add to the assumptions in theorem 18 
the assumption that Stab(ipx) is unitary. Denote a p = 
pi —pG {92), where Q{-) is the Stab(ip) -twirling operation 
(see Eq. (|)J. Th en. 

Pmax(0i —> ^2) = max |p o~p separable^ . 

Remark. Note that a v = Q (pi — pp-z) since we assume 
that Stab(ipi) is unitary. Moreover, if p is small enough 
(see lemma nj) then p\ — pp2 is a positive semi-definite 
separable matrix, and therefore Q {p\ — pp-i) is also pos- 
itive semi-definite separable matrix. This observation 
provides a lower bound on P max , which will be discussed 
after the proof. 

Proof. Based on theorem [18] we require that 
A^ 1 _ > ^, 2 ({pk}) is separable operator. However, since 
G(pi) — Pi (see proposition 11 ) we get that 



= Pi ' 



({Pk})) = G(pi)~pG (p 2 ) 
pG (p 2 ) = Avi->0 2 ({p/™})) 



(19) 



is also a separable matrix. Here, p = 53fc=i Pk ano - {p/ m } 
is a set of m probabilities all equal to p/m. This com- 
pletes the proof. □ 



3. Lower and Upper bounds 

The formula in theorem [18] for the maximum prob- 
ability to convert one state to another is very hard to 
calculate especially if the stabilizer of 0i is not unitary. 
One of the reasons for that is that the optimization in 
Eq. (17) needs to be taken over all sets {pk} such that 
({Pk}) is separable; but to check whether a mul- 



A 



tipartite mixed state is separable or not is an NP hard 
problem [30j . Therefore, it would be very useful to find 
lower and upper bounds to P max in terms of much sim- 
pler expressions. Fortunately, we were able to find such 
expressions. 



An Upper Bound 



Corollary 22. The maximum probability to convert ip\ 
to "02 is bound from above by: 



-Pmax (01 -> 02 ) < 



■MPi)] 



where A max [-] is the maximum eigenvalue of the matrix 
in the square parenthesis. 

Note that another upper bound can be found by re- 

({Pk}) 



placing the requirement in theorem 18 that A 



is separable, with the weaker requirement that it is only 
positive semi-definite matrix. However, this will still lead 
to a non-trivial optimization problem. 



A Lower Bound 

The maximum probability to convert one state to an- 
other is always greater than zero if both states belong 
to the same orbit. In the following corollary we find a 
positive lower bound for this probability. 

Corollary 23. Notations as in theorem \18\ The maxi- 
mum probability to convert 0i to 02 is bound from below 
by: 



c (0! -> 2 ) > 



1 



where A max [-] is the maximum eigenvalue of the matrix 
in the square parenthesis. 

Remark. Since the ADM is defined up to a conjugation 
by an element in the stabilizer group, we can replace p 2 
in the equation above by U^piUl for any Uk G Stab(ip). 
Since A max [pT/ 1 /^] is in general not invariant under such 
a conjugation, we can improve the bound by taking an 
optimiztion over all m elements of the stabilizer; that is, 



Pmax (01 -> ^2 ) > max 
k— 1.....T 



Pl'ulpiUk 



Note also that 



A 



> 



Amax [pi] 



max [Pi P2| 



[P2\ 



with equality if [pi,P2[ = 0. Hence, the maximum eigen- 
value of the associate density matrix can only increase 
under SEP. This implies, for example, that — log A max [pi] 
is an entanglement monotone. In the bipartite case it is 
called entanglement of teleportation [3T] . 



Notations as in theorem |18| A simple upper bound to 
P m follows immediately from corollary |19| The require- 
ment that the matrix in Eq. ( 18 ) is separable implies that 
the matrix is also positive semi-definite. We therefore 
have the following upper-bound: 



D. Infinite Stabilizer Group 

In this section we generalize some of the results ob- 
tained in the previous sections for the case where the sta- 
bilizer group is not finite. Although, for most states the 
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stabilizer group is finite, for many important states, like 
the GHZ type states, W type states, or Cluster states, the 
stabilizer group has non-zero dimension. Furthermore, as 
we will in the examples of the next section, the larger the 
stabilizer group is, the more conversions between multi- 
partite states can be achieved by local separable oper- 
ations. In this context, the fact, for example, that the 
maximally bipartite entangled state can be converted by 
LOCC to any other bipartite state is a consequence of a 
very large stabilizer group (see Eq. ([6])). 

When the stabilizer group of a state ip is infinite, it 
is impossible in general to find a state in the orbit of 
tjj for which the stabilizer group is a unitary group. A 
simple example for that can be found in the stabilizer 
of the three qubits GHZ state discussed in the next sec- 
tion. Nevertheless, theorem 12 in section V A can be 



easily generalized to infinite groups, but the conditions 
are more complicated. 

Theorem 24. Let ) = 9i IV'VHsiV'll an d \ip2) — 
92 1 i>) I II 92ip II be t wo entangled states in the same orbit 
OtU. Then, \ipi) — ¥ |-0 2 ) by SEP if and only if there exists 
to G N and a set of probabilities Po,Pi, ■ ■■,Pm (i-e- Pk > 
and J2T=iPk — lj such that 



'^2pkSlp 2 S k = pi 



(20) 



fe=l 



for some Sk £ Stab(tp), where pi = p^(ipi) for i = 1,2. 

The proof of the theorem above follows the exact same 
lines as in the proof of theorem p~2| Similarly, theorem |l5| 
can also be generalized to non-compact groups in this 
way. However, since the set of to elements, {Sk}, is only 
a finite subset of the infinite group Stab(tfj), it can be very 
difficult to find such a finite subset of Stab(ip) that sat- 
isfies Eq.([20|. Moreover, note that according to propo- 
sition [7] the Stab(ip) is also non-compact which makes it 
impossible to generalize the stabilizer twirling operation 
(see Eq.Q) to this case. Nevertheless, we will see in the 
3 qubits example of the next section that the twirling op- 
eration in Eq. |5| can be a very useful tool in determining 
local conversions of multipartite states with non-compact 
stabilizers. The idea to define the group T as the inter- 
section of Stab(ip) with the group of local unitary matri- 
ces is useful in general (not only for the 3 qubits example 
below), but we leave the general analysis for future work. 



VI. EXAMPLES 

In this section we give examples demonstrating the ap- 
plications of the main results to specific cases. We start 
by showing how our techniques produce all the known 
results about the necessary and sufficient conditions for 
local interconversion among pure bipartite states. Then, 
we focus on Hilbcrt spaces of three and four qubits. We 
start with the four qubits case because the stabilizer of 
generic states in four qubits is finite, while the stabilizer 



of the three qubits GHZ state turns out to be a two di- 
mensional non-compact group, and therefore the analysis 
is a bit more complicated. 



A. The Bipartite Case 

In this subsection we consider the bipartite example 
with a Hilbert space H2 = C d <E)C d . Up to local unitaries, 
the only critical state is the maximally entangled state 

d-i 



£iw>, 



i=0 

where is an orthonormal basis in C d . Clearly, in 

this case the orbit is dense and consists of all states 
in C ®C with Schmidt number d. In order to apply the 
theorems from the previous sections, we need to calculate 
first the stabilizer group. 

Lemma 25. The stabilizer group of ip is given by 
Stab{ip) = {s- 1 ®S T Se GL(C d )} , 

where S T is the transpose matrix of S . 

Proof. It is straightforward to check that ® S T \ip) = 
\tp). Now, if A ® B e Stab(tp) then 

= A®B\xp) = (A®B) (A- 1 ® A T ) = I®BA T \ii) . 

Thus, BA T = I which implies that A ® B is of the form 
given in the lemma. □ 

Clearly, the lemma above implies that Stab(ip) is not 
compact. However, in |32j it was shown that when con- 
sidering only pure bipartite states, any LOCC transfor- 
mation can be achieved by the following protocol: Alice 
performs a quantum measurement and send the result to 
Bob. Based on this result, Bob performs a unitary ma- 
trix. Hence, all LOCC protocols can be achieved by the 
following completely positive map. For all \ijj) € C d (S>C d , 



A (10X01) 



52Mk\4>)(Wl 

k 



M k = A k 



Bi 



where B^ is a unitary matrix. Thus, without loss of 
generality, for the bipartite case we can replace Stab(ip) 
in theorem 24 with the subgroup T C Stab(ip) given by 



V 1 



Ve 



U(d)} 



Note that T is a compact unitary group. 

Now, let and \1fj2) be as in theorem 12 applied to 



the bipartite case, and denote by p\ and pi their ADMs, 
respectively. Towards the end of Sec. |IV| we have seen 
that the ADMs p\ and p 2 can be written as 



Pi 



Pi l) 



for i = 1.2 
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where p^ 1 and p)?' are the the reduced density matrices 
of |V>i) and |t/> 2 ), 
that the map |^x! 
if and only if their exists probabilities {pk} and elements 
5 fe G T such that 



respectively. Thus, theorem 24 implies 
■ — > \ip2) can be implemented by LOCC 



^^PkSlp 2 Sk = pi ■ 

k=l 

But this condition is equivalent to 

m 



k=l 



where {Vfe} is a set of m unitary matrices. This is the 
exact same condition obtained in (27) for pure bipartite 
LOCC transformations and is equivalent to the majoriza- 
tion condition. In the same way, we can show that when 



applying theorem 15 to the bipartite case, we can ob- 
tain the same results found in [28] for non-deterministic 
LOCC transformations. 



B. Four qubits 

In this subsection we consider the Hilbert space of four 
qubits H 4 = C 2 ® C 2 <g> C 2 ® C 2 . The SLOCC group is 
given by G = SL(2, C) ® SL(2, C) ® SL(2, C) ® 5L(2, C). 
The analysis of four qubits SLOCC classes can be found 
in [TO] EU EU |33] (and references therein) , and we refer 
the reader to these references for more details about some 
of the statements we make here. 

In four qubits the critical set can be characterize ele- 
gantly by the following 4 dimensional subspace: 

Crit{H±) = KA, 

where K = SU{2) ® SU{2) ® SU(2) ® SU{2) is the set 
of special local unitary matrices, 



A — ^ZqUq + Z\Ul + z 2 u 2 + z 3 u 3 

with 



z ,z 1 ,z 2 ,z 3 G 



Uq 
U 2 



= \<j> + )\cj> + ) , Ul = \r)\r) 

\1>-)\4>-) 



"3 



and 1^) = (|00>±|ll))/V2and 1^) = (|01>±|10))/V2. 
In [33] the class A was denoted by G abed- 
Note that not all the states in A are generic. For exam- 
ple, the state uo is not generic. It is possible to show El] 
that a four-qubit state <j> is generic if and only if (a) the 
dimension of the (non-normalized) orbit G(f> is maximal 
(i.e. 12), and (b) the orbit Gcj> is closed in Ha- 

The set of all generic states in Crit(Hi) is given by 
KT>, where 



V 



ZqUq + Zlttl + Z 2 U 2 + Z3U3 



^ z) for i^i} 



In the following theorem we prove that for all these states 
the intersection of stabilizer group with G is the Klein 
group, consisting of four elements. 

Proposition 26. Let ip G T> and denote by St(ip) the 
intersection of Stab(tp) with G (see corollary\10\j. Then, 



St(1>) = {I,X,Y,Z}, 

where X , Y and Z are the three 2x2 Pauli matrices and 
C = C ® C ® C <g> C for C = X,Y, Z . 



, 11 " ,l,L 7^ z j f° r 



Proof. If <f> G A and cf> — Ylj=o z j u j an d z, 
i 7^ j then we say that </> is generic. If g G G and 
A then if /3 = ^ x^Ui then g/3 — ^ EiX a iUi with cr a 
permutation, £j = ±1 and £o£i£2£3 = 1- We will show 
that the stabilizer in G of a generic element in A has the 
property that it leaves A invariant. We therefore see that 
if g<f> = cj) with generic then g\A = I. We now give more 
details (see [2"5j). 

To complete the discussion we reformulate the action 
of G on Ha- We first consider SX(2,C) <g> SX(2,C) act- 
ing on %2- There is a non-degenerate symmetric form, 
(..., ...) on n 2 that is invariant under SL(2,C)®SL(2,C) 
(take the tensor product of the symplectic form on C 2 ). 
Relative to an orthonormal basis of T~L 2 the image of 
SL(2,C)®SL(2,C) is 50(4, C). Thus the action of G on 
H4, ='H 2 ® H 2 is equivalent to the action of S0(4, C) <X> 
50(4, C) on C 4 ® C 4 . We can use the form (..., ...) to 
identify C 4 ® C 4 with M 4 (C) with 50(4, C) x 50(4, C) 
acting by (g,h)X = gXh T — gXh" 1 . One can check 
that with these identifications we can take A to be the 
diagonal elements. If we consider the group 50(8, C) 
then we can consider G to be the subgroup 



9 
h 



| 5 ,ftG50(4,C) 



as 



The space H4 can be imbedded in Lie(50(4,( 

X T 'I A' G .V„ 



The space A now corresponds to a Cartan subalgebra of 
Lie(SO(8, C)) and the assertion above is implied by the 
fact that the Weyl group of 50(8,C) is given as above. 
This implies that the stabilizer of a regular element in A 
fixes every element in A. In particular if (g, h) is in the 
stabilizer then 

gIh T = I 

so g = ft. Also if for all X we have gXg^ 1 = X then 
g is diagonal. Since g is orthogonal gg T = I so g 2 = I. 
Hence g is diagonal with entries ±1. Now unwind the 
identifications. □ 

In 4-qubits there are 4 generating SL invariant poly- 
nomials [10]. One of these generating polynomials is of 
order 2 and is given by: 

f(l/j) = (lp*\0-y <g) a y <E> (Ty <X> 0-y\lp) . 
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The absolute value of this SL-invariant polynomial is 
called the 4-tangle. If if) € Criti^i^) then 



Z Q + z l + z 2 + ^3 



Now, from proposition [8] (see also the remark after the 
proof of proposition [8]) it follows that if f(ip) ^ then 

Stty) = Stab(ip) . 

That is, for all the generic states in V with f(ip) ^ 0, the 
stabilizer group is the 4 element Kle in gr oup Stabfy) = 
{I, X ,Y , Z}. Hence, the results in VA implies that a 



state ip e T> with f(tp) ^ 0, can be converted to 
\g\ip) by SEP, if and only if 



9*9 



9{?)=I 



where Q is the Sia6(^)-twirling operation defined in 
Def.[3] Moreover, it is possible to characterize completely 
all the ADMs a that satisfy the equation above. This 
can be done by applying the factorization in Eq.Q to 
the Klein group. 

For states ip £ T> with f(ip) — one can use corol- 
lary |10| to find out the stabilizer group. For example, the 
state (see [TU] for important properties of this state) 

\L) = -j= (u + ujux + uu 2 ) , 

is in D and satisfies f{\L)) = 0. In fact it is possible to 
show that it vanish on all ST-invariant polynomials of 
degree less than 6, and that there exists an SL- invariant 
polynomial of degree 6 on which it does not vanish [TU] , 
Hence, from corollary [Tojit follows that Stab(\L) ) / St(\L) ) 
is a cyclic group with three elements. The generator of 
the group is given by g = ujp 3 , where p 3 € K is the local 
unitary permutation that takes uq — !• u 2 , u\ —t uq, U2 — > 
u\ and M3 — > U3. Hence, now that we found Stab(\L)), 
we can apply all the theorems in the previous section to 
determine the states to which \L) can be converted by 
SEP. 



C. Three Qubits 

In this subsection we apply the results from the previ- 
ous section to the case of three qubits. In particular, we 
will consider local transformations between two states in 
the GHZ SLOCC class [9]. Since the stabilizer group of 
the GHZ state is non-compact, we will only be able to 
solve the problem partially. 

We denote the GHZ state by <\> = ^(|000) +|1H)). A 
straightforward calculation shows that the stabilizer of <j> 
in G is the group 



st{4>) 





,-1 



, S1S2S3 — 1 



Note that if t = t± ® to, ® t$ as above then we have 
i|000) = |000) ,t |111) = I 111) 

t|001) = S3- 2 |001) ,t |010) = S2 2 1010) ,t 1 100) = s^ 2 l 10 °) , 
* |011) = si |011> , * ]101) = s 2 2 1 101) , i 1 110) = si |110> . 

Now, the 3-tangle of the GHZ state is non-zero. Since 
the tangle is the absolute value of an homogeneous SL- 
invariant polynomial of degree 4, from corollary [TU] it 
follows that Stab((f>) / St(tfi) is a cyclic group with two 
elements. The two elements are 

Stab{(j>)/St{4>) ={I , X ® X ® X} , 

where X is the x-Pauli matrix. Stab(<f)) is therefore the 
union of St(<j>) and hSt{<j)), where h = X ® X ® X. 

Set T equal to the intersection of Stab(<fi) with the 
group of local unitary transformations, and To the inter- 
section of St(4>) with the group of local unitary transfor- 
mations . The T-twirling operation is therefore 



7» 



-h%{a)h 



where 



To(o-) = J dt t f at 



with dt — -^jdQxdQi for Sj 



6, = 0. 



Now, from theorem 24 it follows that if tpi and ip2 are 
two states in the GHZ class, then it is possible to convert 
-01 to ip2 by SEP if and only if there exists probabilities 
{pk} and elements € Stab{4>) such that 



^2Pk9kP29k 

k 



Pi 



where p\ and p 2 are the ADMs of i\>\ and ^2 , respec- 
tively. In order to get a simpler condition we apply the 
To-twirling operation on both sides of this equation. This 
gives, 



E 

k 



Vk%(g\p29k) = % (pi) . 



Since St(cf>) is commutative and hSt((f>) = St(<j))h we get 
]pkgt%{p2)9k = T Q (pi) . 



k 



Now, without loss of generality assume gk £ St(<j)) for 
k < a, and gk = hfk with fk G St(<fi) for k > a. Then, 

^2pkgtT (p2)gk = ^2pkgtgkTo(p2)+^2pkftfkT (hp 2 h) 



k<a 



k>a 



where we have used the fact that To(p 2 ) commutes with 
every element of T , and since St(<p) is the complexi- 
fication of To we see that 7o(p2) commutes with every 
element of St(tp). This implies that 



k<a 



Pkg k gkTo(p 2 ) + y2pkf k fkT (hp 2 h) = % (p x ) 
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In the case where 7o (pi) — pi, the above condition is 
both necessary and sufficient. In particular, if we take 
tpi to be the GHZ state <f>, then p\= I and we get 



k<a 



Pkgtgk%{p 2 ) + y2pkftfkT (hp 2 h) 



k>a 



The above condition provides necessary and sufficient 
condition for local separable conversion of the GHZ state 
ipi = 4> to some other state ^2- In particular, the condi- 
tion implies that the GHZ state can be converted by SEP 
to any state ip2 with T(p2) = I since we can take in this 
case and fk to be unitaries. However, it is left open if 
there are other states that can be obtained by SEP from 
the GHZ state. 



VII. CONCLUSIONS 

Quantum information can be viewed as a theory of 
interconversions among different resources. Multipar- 
tite entangled states, like cluster states, are resources for 
quantum information and quantum computing. In this 
paper we considered transformations that can be imple- 
mented on pure multipartite entangled states by local 
separable operations. Quite remarkably, we were able 
to generalize both Nielsen majorization theorem [2 7) and 
Jonathan and Plenio theorem [28] for transformations in- 
volving multipartite pure states. This generalization was 
expressed in terms of the the stabilizer group and the as- 
sociate density matrices (ADMs) of the states involved. 

For the bipartite case with a Hilbert space C d (g> C d , 
both the majorization theorem in |27| and the conditions 
in [55J can be expressed in terms of Vidal's entanglement 
monotones [35] defined by 



E k (i>) = J2 X k k = l,2,...,d 



where Afc are the ordered (from the biggest to the small- 
est) eigenvalues of the reduced density matrix, and ip is 
a bipartite state. Here we were able to extend this def- 
inition to the multipartite case by taking the AfeS to be 
the ordered eigenvalues of the ADM of ip- However, un- 
like the bipartite case, for local transformations involving 
multipartite states, the conditions given in terms of these 
functions are only necessary but not sufficient. It is there- 
fore impossible to express the necessary and sufficient 
conditions only in terms of functions of the eigenvalues 
of the ADMs involved in the transformation. Instead, we 
express the conditions for local transformations in terms 
of matrix equalities given in theorems [12] and |15[ 



The necessary and sufficient conditions gets extremely 
elegant for transformations that involve the maximally 
entangled state of the SLOCC orbit. That is, every 
SLOCC orbit with a finite stabilizer group has a max- 
imally entangled state (this is usually a critical state, i.e. 
a normal form) . Its stabilizer group consisting of unitary 
matrices, and the necessary and sufficient conditions to 
convert this state to another state in its SLOCC orbit, 



are given in Eq.(10l and expressed elegantly in terms of 
the stabilizer twirling operation (see definition [3b. This 
simple expression made it possible to characterize all the 
states to which a maximally entangled state can be con- 
verted by SEP. 

Both theorem [12] and [15] assume that the stabilizer 
group of the state representing the SLOCC orbit is fi- 
nite. We also argued that for a dense set of states in 
H n with n > 3 the stabilizer is finite. Therefore, these 
theorems hold true for most states in Tl n . Moreover, 
from continuity arguments, it follows that the necessary 
conditions in Eq.(13) given in terms of the continuous 
functions Ek(ip) hold true for all states in T-L n . That is, 
the necessary conditions that generalize the majorization 
theorems in [57] and [55] hold true for all the states in 
the Hilbert space. In particular, this implies that Ek are 
entanglement monotones. 

When local deterministic transformations are not pos- 
sible, we found an expression for the maximum probabil- 
ity to convert one multipartite state to another by SEP. 
This expression can be simplified dramatically if one of 
the two states involved is critical. However, since in gen- 
eral the expression is relatively complicated, we found 
simple lower and upper bounds for the probability of lo- 
cal conversion. 

There are many important SLOCC orbits, like the 
GHZ class, that have a non-finite (in fact, non-compact) 
stabilizer group. For such orbits we also found necessary 
and sufficient conditions for a local conversion among 
states in the same orbit, but these conditions are more 
complicated since the stabilizer group is not unitary. We 
are willing to conjecture that for such orbits theorems [57] 
and [28] hold true if one replace the non-compact stabi- 
lizer group, with the subgroup obtained by the intersec- 
tion of the stabilizer with the local unitary group. 

There are many other interesting related problems that 
where not discussed here. One such problem concerns 
with the asymptotic rates to interconvert one multipar- 
tite state to another. We hope the work presented here 
will be useful for such a research direction. 
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